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We consider a three-dimensional lattice (7(1) x (7(1) superconductor in the London limit, with 
two individually conserved condensates. The problem, generically, has two types of intercomponent 
interactions of different characters. First, the condensates are interacting via a minimal coupling to 
the same fluctuating gauge field. A second type of coupling is the direct dissipationless drag repre- 
sented by a local intercomponent current-current coupling term in the free energy functional. The 
interplay between these two types of interactions produces a number of physical effects not present 
in previously investigated U{1) x U{1) models with only one kind of intercomponent interaction. 
In this work, we present a study of the phase diagram of a f/(l) x 17(1) superconductor which 
includes both of these interactions. We study phase transitions and two types of competing paired 
phases which occur in this general model: (i) a metallic superfluid phase (where there is order only 
in the gauge invariant phase difi'erence of the order parameters), (ii) a composite superconducting 
phase where there is order in the phase sum of the order parameters which has many properties 
of a single-component superconductor but with a doubled value of electric charge. We investigate 
the phase diagram with particular focus on what we call "preemptive phase transitions" . These 
are phase transitions unique to multicomponent condensates with competing topological objects. A 
sudden proliferation of one kind of topological defects may come about due to a fiuctuating back- 
ground of topological defects in other sectors of the theory. For (7(1) x (7(1) theory with unequal 
bare stiffnesses where components are coupled by a non-compact gauge field only, we study how 
this scenario leads to a merger of two t/(l) transitions into a single (7(1) x (7(1) discontinuous 
phase transition. We also report a general form of vortex-vortex bare interaction potential for an 
N-component London superconductor with individually conserved condensates in the presence of 
interactions mediated by a fluctuating gauge field as well as intercomponent drag. 

PACS numbers: 67.85.De,67.85.Fg,67.90.-|-z,74.20.De,74.25.Uv 



I. INTRODUCTION 

Phase diagrams and critical phenomena in superfluids 
and superconductors with U{1) symmetry are well un- 
derstood theoretically and well investigated numerically. 
The understanding is largely based on identifying and de- 
scribing the behavior of proliferating topological defects. 
In two dimensions, a transition from a superfluid to a 
normal state can be described as unbinding of vortex- 
antivortex pairs, which disorders the superfluid phase 
yielding a Berezinskii-Kosterlitz-Thouless transition into 
a normal state. ^ In three dimensions, the topological de- 
fects of (7(1) theory are vortex loops, proliferation of 
which yields a continuous phase transition in the 3Dxy 
universality class in the case of superfluids (with global 
(7(1) symmetry), or inverted 3Dxy in the case of super- 
conductors (with local U{1) symmetry) .^^Hil However, it 
was recently found that in interacting mixtures of (7(1) 
symmetric condensates the situation changes principally, 
yielding much more complex physics, different phase dia- 
grams and transitions. Many aspects of the phase transi- 
tions in systems with several interacting components are 
still poorly understood and debated. 

The main important new aspect arising in an inter- 
acting mixture is connected with the fact that, as re- 
viewed below, under certain quite generic conditions the 
vortices with high topological charge (or bound states of 
vortices) acquire crucial importance for various aspects in 



the physics of these systems. This is in contrast to single- 
component systems where only the lowest-topological- 
charge defects (i.e. only vortices with 27r phase winding) 
are important. The complexity arising from the rele- 
vance of topological defects with high topological charge 
include formation of what is called "metallic superfluid 
phases", in context of electrically charged systems, or 
"paired phases", in context of electrically neutral sys- 
tems. In these states no conventional real space pairing 
takes places. However, there is order only in the sum or 
difference of the phases of the condensate, with phases 
being individually disordered.'^ltl^ Moreover, it also re- 
sults in a complicated and still poorly understood nature 
of the phase transitions from a fully symmetric state to 
a state with all symmetries broken,'i2Hl3l ^i^n there is a 
competition between proliferating low- and high- order 
topological defects. This is again a phenomenon which 
has no counterpart in single-component systems. Vari- 
ous aspects of related effects were also studied in differ- 
ent models with a compact gauge field and with SU{2) 
symmetry.l^ 

Recently, it has been found that two kinds of inter- 
component interactions lead to the novel mixture-specific 
phenomena mentioned above. Namely, in a mixture of 
charged condensates, the intercomponent interaction is 
represented by the coupling between the charged com- 
plex sc alar mat ter fields mediated by a fluctuating gauge 
figljjIIIHE o | i5 | i6 | ^j^g other hand, in the case of an elec- 
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trically neutral condensate mixture, some related (but at 
the same time principally different) effects can be pro- 
duced by a stron g dissipationless drag (current-current 
interactiorPi^Ii^ which in some physical situations is also 
called Andreev-Bashkin interaction).^* The intercompo- 
ncnt couplings by gauge field and the dissipationless drag 
have so far only been studied separately, while in a 
generic t7(l) x [/(I) system, terms leading to both of these 
effects are allowed by symmetry. Thus, generically the 
phase diagram and critical phenomena in a U{1) x C/(l) 
system is a problem with two coupling constants. The 
interplay between them has, to our knowledge, not been 
investigated so far. 

In this work, we report a quantitative study of a 
generic U{1) x U{1) London superconductor which has 
both kinds of intercomponent coupling (gauge field and 
current-current drag). This includes, in particular, the 
situations where these two different kinds of intercompo- 
nent couplings compete with each other. 

This paper is organized as follows. In Sec. [llj the 
general model we consider is introduced, and the neutral 
and charged modes and the vortex representation of the 
general model, obtained by a duality transformation, are 
identified. Sec. IIIII is devoted to the numerical methods 
we employ in this study. The results obtained in the 
special case with no intercomponent dissipationless drag, 
is presented in Sec. |IV[ followed by the results of the 
general model with competing gauge field and Andreev- 
Bashkin interactions in Sec. |V) The conclusions are given 
in Sec. |VI[ We also present analytical details presenting 
the duality transform for a general A'^-component model 
in Appendix [A| and a derivation of the expression for the 
gauge field correlator in Appendix [B} 



II. THE MODEL 

We study a generic two-component London supercon- 
ductor. In the London limit, one neglects the fluctua- 
tions of the density fields \ipi.2\ of the complex scalar 
functions ipj = \ipj\e^^^ describing two superconducting 
components (i.e. setting |V'i,2| ~ const). Fluctuations of 
the phases 0j, and the gauge field A are allowed. The 
compact support of the phase-variables 9i € [0, 27r) im- 
plies that phase-fluctuations lead to vortex-excitations, 
capable of destroying superconductivity/superfluidity, in 
this system. The London limit is an adequate approxima- 
tion for many properties of strongly type-II superconduc- 
tors, and in fact transcends the validity of the Ginzburg- 
Landau theory. The free energy density of this system 
can be written as 



Pd 



(1) 



{V0i - eiA- V6I2 + e2A)^ 



between the two charged condensates via a fluctuating 
gauge field A, we include a direct intercomponent dissi- 
pationless current-current interaction with strength p^, 
which has the forrrfl^ 



Fdrag = Pd^^Oi ~ CiA) • {^02 - GaA). 



(2) 



It is a part of the last term in ([T]). The particle currents 
of both species then depend on the common vector po- 
tential and superfluid velocities of both condensates (i.e. 
particles belonging to one condensate can be carried by 
superfluid velocity of the other), 

ji = (Pi - Pd)(V0i - eiA) + pa{V92 ~ e2A), (3) 
h = (P2 - Pd){^02 - 62 A) -f pdiy0i - eiA). (4) 

For generality, we allow for unequal charges in the 
two-components of the system, examples of the systems 
with oppositely charged condensates are given below. 
Note that the drag term implies that there is a stabil- 
ity criterion that must be applied to the system. If 
exceeds a critical limit, to be determined below, the spec- 
trum of the system will be unbounded from below and 
hence the theory will be ill-defined. The bare stiffness co- 
efficients Pj must be positive, pj > 0, on simple physical 
grounds. 

The physical model in ([T]) is discussed in the context of 
the projected quantum ordered states of hydrogen or its 
isotopes at high compression^"^ •^Q- where the differ- 
ent fields correspond to condensates formed by electrons, 
protons or deuterons. A similar model appears in some 
models of neutron stars interior where the two fields rep- 
resent protonic and S~ hyperon Cooper pairs.^^ More- 
over, the model with equal phase stiffnesses pi = p2 and 
charges ei =62, appears as an effective model inthe the- 
ories of easy-plane quantum antiferromagnets.'22l2ll Re- 
lated mode ls we re also studied in various contexts in low 
dimensions 

The model has topological excitations which are vor- 
tices with 27rnj, rij = ±I,±2, ... phase winding in the 
phase of component j. We denote vortices by the pair 
of integers (77,1,77.2) characterizing phase windings of the 
vortex in question. Thus, vortices with phase winding 
in only one component are denoted (1,0) or (0, 1). The 
model also possesses composite vortices where both inte- 
gers associated with the phase windings (around or nearly 
around the same core) in the two species are nonzero. In 
this paper, we will only consider the composite vortices 
(1,1) and (1,-1) which have co-directed and counter- 
directed phase windings in the two components, respec- 
tively. However, composite vortices with higher topolog- 
ical charges, such as ( 1 , 712 ) o r (tt-i,!), may be relevant 
under certain conditions .'iSESl 



A. Charged and neutral modes 



where pj physically represent the bare phase stiffnesses of 
the problem. In addition to the intercomponent coupling 



By separation of variables ,1^122121] .^g rewrite the 
model in ([T]) in a form where the composite charged and 
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neutral modes are explicitly identified, 



B. 



The case pd = 0, ei — €2 — e 



F = 



P1P2 - PdiPl + P2) 



(62^01 - eiV6'2)^ 



2 A ^2 



(V X A)^ 



(5) 



where the coefBcients are given by 



Pi - Pd 
Pd 



Pd 

P2 - Pd 



(6) 



and 



[pi- Pd)el + {p2- Pd)el + 2pdeie2. (7) 



Throughout the paper, there is an implicit sum over re- 
peated component indices. The coefficient mo should not 
be confused with the mass of the components. These are 
included in pj, whereas mo determines the inverse bare 
screening length of the screened interactions in the sys- 
tem, details will be given in Sec. |II C[ The first term of 
([5]) is identified as the neutral mode that does not couple 
to the vector potential. The second term is the charged 
mode, characterized by its coupling to the vector poten- 
tial. 

From Eq. it is seen that for stability of the sys- 
tem (in the sense that the free energy functional should 
be bounded from below) the coefficient of the first term 
should be positive. It is readily shown that the criterion 
for this is that 



Pd < 



Pi P2 
Pi + P2 



(8) 



Note that this criterion is identical to the one derived in 
Ref. [131 and does not depend on charge. Actually, there 
are no restrictions on the value of the electric charge e, 
to obtain a well-defined theory. 

Note that in ([T]), the phases of the two components 
do not represent gauge invariant quantities. However, 
when the model is rewritten on the form in ([5]), observe 
that the neutral mode identifies a linear combination of 
the phase gradients that is a gauge invariant quantity 
decoupled from the vector potential A, V(e2(?i — ei02)- 
Thus, the U{1) x U{\) symmetry of the model may be 
interpreted as possessing a "composite" electrically neu- 
tral (or "global") U{1) symmetry associated with the 
phase combination of the neutral mode, and a "com- 
posite" U{1) gauge symmetry which is coupled to vec- 
tor potential A and thus is associated with the charged 
mode. Importantly, the identification of the charged and 
a neutral mode does not imply that the modes are de- 
coupled, because both modes depend on phases 0i which 
are constrained to have 27r x integer phase windings. 



We briefiy review the physics of a two-component 
C/(l) x U(\) superconductor with individually conserved 
condensates, coupled only by the gauge field i.e. in the 
absence of Andreev-Bashkin (i.e. mixed gradient) terms. 
In the London limit the free energy may be read off from 
Eq. ([5|, 



1 



P1P2 



fV( 



2pi+ P2' 
1 [piV^i 



P2V02 - e(/9i +P2)A]' 



(9) 



Pl+ P2 



^(Vx A) 



The important new physics arising in the model ^ 
compared to single-component GL model is that the 
lowest-order topological defects with a 27r phase wind- 
ing only in one phase 9i have a logarithmically diverging 
energy per unit length, while vortices where both phases 
have 2tt winding have finite energy per unit length j^'^^^ 
Under certain conditions vortices where both phases 
wind, i.e. (1,1), can proliferate without triggering a pro- 
liferation of the simplest vortices (1,0) and (0,1). 

Consider now a composite (1,1) vortex. Such an exci- 
tation, if vortices in two components share the same core, 
has nontrivial contribution to the following terms in the 
free energy ^ 



poff 

-^(1,1) 



1 [pi V6I1 + P2W02 - e{pi + p2)Af 



Pi + P2 



1 



(10) 



(V X A) 



If the (1,1) vortex has phase windings around a common 
core, it can be mapped onto a vortex in a single compo- 
nent superconductor. Then, by increasing electric charge 
one can make the energy cost of a vortex per unit length 
in a lattice London superconductor arbitrarily small (be- 
cause the vortex energy depends logarithmically on the 
penetration depth which is in turn a function of electric 
charge). Thus, in a lattice London superconductor the 
critical temperature of proliferation of the vortices can be 
arbitrary small if the value of the electric charge is suf- 
ficiently large. Therefore, in the two-component model 
([9]) one may, by increasing the value of electric charge, 
proliferate (1,1) vortices without proliferating individual 
vortices (1,0) or (0,1). The latter two produce a phase 
gradient in the gauge invariant phase difference 0i — 02- 
This features a stiffness which is not renormalized by the 
proliferation of the (1,1) vortices. 

Since the (1,1) vortices do not have a topological charge 
in the phase difference, they cannot disorder the first 
term in but they disorder the charged sector repre- 
sented by the second term. The resulting state therefore 
features long-range ordering in the phase difference and 



can be characterized by (e 



ii0l 



2)\ 



^ 0, while (e 



0, 
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(jgi6>2^ = 0, and there is no Meissner effect. The free en- 
ergy for the resulting phase is given by the following term 
(i.e. it has only broken global U{1) symmetry) while the 
stiffness of the charged U{1) mode is renormalized to zero 
by proliferated composite vortices, 



f: 



lOff 



1 



PlP2 



2 Pi + P2 



[V(01 - 02)? 



(11) 



The proliferation of composite defects resulting into this 
state was shown to arise in two-dimensional systems at 
any finite temperatures.^^In three dimensions, this phase 
can be induced by a mag netic field via melting of a 
composite vortex lattice.'^ An analogous phase was also 
found in a three dimensional lattice superconductor aris- 
ing without applied external field from fluctuations if the 
value of the electric charge is very large.ESl Since there is 
no Meissner effect in the resulting phase, but at the same 
time there is a broken neutral U{1) symmetry, the term 
metallic superfiuid (MSF) was coined for it-^^Also related 
phases are sometimes called "paired phases" ^ The lat- 
ter term is motivated by the fact that in such situations 
the (quasi-) long-range order is retained only in some 
linear combination of phases while individual phases are 
disordered. Importantly it should not be confused with 
the conventional "real-space" pairing of bosons. 



1. The case pi = p2 

Consider the case where pi — P2- At high values of the 
electric charge e, the model was shown to feature a MSF 
phase without applied field.'^ This implies that at large 
e the system undergoes two phase transitions when the 
temperature is increased. The first transition is from a 
state with broken U{1) x U(l) symmetry into the MSF 
with broken U{1) symmetry, driven by a proliferation of 
composite (1,1) vortices. The second transition is one 
where the remaining broken U{\) symmetry is restored 
by proliferation of individual vortices, resulting in a nor- 
mal state. At low values of e, one cannot separate char- 
acteristic temperatures of the proliferation of composite 
and individual vortices and thus, the model should have 
only one phase transition from broken U{1) x U{\) to a 
normal state. In the case pi — P2, the latter phase tran- 
sition was conjectured to be a continuous phase transi- 
tion in a novel universality class in the work of Ref. fT^ 
However, subseq uent works show that the phase transi- 
tion is first order .'iSmi^ see also Moreover, the analysis 
performed in Ref. ITUl indicates that the U{1) x U{1) 
to a normal state transition is first order for any val- 
ues of electric charge in the pi — p2 model. Note that 
the standard theories o f vo rtex loop proliferation yield a 
second order transition.!^ An analysis of a simpler two- 
component model (with no gauge field coupling, but with 
direct current-current coupling) which, like the model 
([9]), also features low-energy composite vortices, provides 
some evidence that the first order transition takes place 



whenever a restoration of the U{1) x U{1) broken sym- 
metry is driven by proliferation of competing tangles of 
different kinds of vortices^ e.g. tangles of (1,0), (0,1) 
vortices and a tangle of (1,1) vortices. The term "pre- 
emptive vortex-loop proliferation transition" was coined 
for this scenario .'^ Note that in a charged U{1) x U{1) 
theory for arbitrary values of electric charge one cannot 
rule out in a simple way that composite vortices partici- 
pate in a competition with the individual vortices in the 
symmetry-restoration transition since composite vortices 
have finite energy per unit length. 



C. Dual model 

We will now perform a duality transformation that re- 
duces the model in ([T]) to a theory of interacting vor- 
tex loops of two species. These are the topological ob- 
jects which drive the phase transition between the normal 
state and a state with broken symmetries in the systems 
we consider. When the phases and gauge field are fiuctu- 
ating the statistical sum of the London two-component 
superconductor with intercomponent drag can be repre- 
sented as follows 

Z = J T^^i J ^^^2 J VA e"'^, 

{[VxA(r)]2 



S 



d-'r 



+ [V0,(r) - e,A(r)] R,k [V0fc(r) - efcA(r)] }, 



(12) 



where /3 is the inverse temperature. 

We now choose the gauge V • A(r) = and Fourier 
transform the action. The action is then written as 



j(q) Rjk o , '2 ) Ufc(-q) 



(13) 



where the Fourier transform of W9j{r) is denoted by 
Uj (q) . Moreover we have completed the squares of the 

gauge field with A(q) = A(q) - ejRjk'Uk{q)/{q'^ + ml) 
as the shifted gauge field. By integration of the shifted 
gauge field, the model is written 



d^gU,(q) i?, 



(14) 



with the phases as the only remaining fluctuating quan- 
tities. The phase gradient can be decomposed into a 
longitudinal and a transverse part, Uj(q) = [Uj(q)]L -I- 
[Uj(q)]2^, where the longitudinal component corresponds 
to regular smooth phase-fluctuations with zero curl, i.e. 
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"spin waves". Hence, the longitudinal part is curl- free, 
q X [Uj(q)]L = and the transverse part is divergence- 
free, q- [Uj(q)]r — and thus it is associated with quan- 
tized vortices. One can introduce the field mj(q) which 
is the Fourier transform of the integer valued vortex field 
for component j. 



iq X [Uj(q)]T = 27rmj(q), 



J 



1, 2. 



(15) 



Note that this relation yields the constraint q-mj(q) = 0, 
i.e. the thermal vortex excitations in the theory are 
closed loops as required by the single- valuedness of the or- 
der parameter in an infinite system. In the following, we 
will disregard the longitudinal phase-fluctuations since 
the physics at the critical points in this system is gov- 
erned by the vortex excitations and not the spin-waves. 
The latter are known to be innocuous and incapable 



of destroying long-range order in three-dimensional sys- 
tems. By (15 1, the transverse phase gradient is explicitly 
written 



[lJM)]T-2m^^^'^, J = 1,2, (16) 



and thus, we finally express the statistical sum via vortex 
fields 

mi m2 

S = 2l3Tr^ J d^q m, (q)F,fe(g2)mfe(-q). (17) 
The vortex- vortex interactions are given by 



I Pi- Pd- \A 



VA^) - ( R,. - ^J^H^M^^] ^ V{q^) 



{sjRji 



Pd 



ejefc-Rji-Rfc2 ^j^k^jiRk^ \ 



ejekRjiRk2 



q'^ + rriQ q^ q^ + itIq 

ej-efc_Rji_Rfc2 ^ ^ (ej-Rj2)^ (ej-Rj2)^ 

P2- Pd ^ 

+ ttIq q"^ q'- + mk ) 



(18) 



Here, we have used the identity 

(2^) 



[U,(q)]T • [U,(-q)]T = ^m,(q) • mfc(-q), (19) 



found by Eq. (16 1. We may now interpret mg, given 
by ([t]), as the inverse bare screening length that sets the 
scale of the Yukawa interactions in the system. 

We remind the reader briefly of what is known for 
the one-component case, i.e. p2 ~ 0, Pd = 0, 62 = 0, 
pi = p ^ 0, ei = e ^ in Eq. (18). Then, we have 
Vii = [p - e2p2/(g,2 TO2)]/g,2^ ^^^Yi ml = pe^. Thus, 
^11 = p/il^ + ''^o) ^ screened interaction between the 
vortices, mediated by the fluctuating gauge field. This 
is drastically diff'erent from the multi-component case, 
where one fluctuating gauge fleld is incapable of fully 
screening interactions between vortex excitations in all 
condensate flelds-I^I^Ell 

The interactions between vortex elements in the sys- 
tem are generally seen to include two parts: A long range 
Coulomb interaction with no intrinsic length scale that 
decays as 1/r, and a short range Yukawa interaction with 
an exponential decay. Note that in the index repre- 



sentation of (18), the first term, Rjk/q" will dominate 



the second term, eiemRijRmk/['f'{<l^ + "^o)] ^ 
short distances when q^ is large, because the Yukawa and 
Coulomb part of the second term will cancel each other. 
Effectively, at short distances, the vortices will interact 
as if the gauge field does not fluctuate. On the other 
hand, at large distances, when q'^ is small, the Coulomb 



part of the second term will dominate its Yukawa coun- 
terpart and the second term will be of the same order as 
the first term ~ q~^ ■ Thus, the l/r contributions from 
the gauge field mediated interactions between vortices 
sets in when intervortex separation becomes larger than 
the characteristic distance TOq ^. Also note that by de- 
creasing the gauge field coupling constant e, ttiq ^ grows 
and so does the distance where the effects of the gauge 
field are negligible. In particular, when = in ( 18 ) 



(this corresponds to the work in Refs. 1519^ . we have the 
case that the interactions between elementary vortices of 
different species tend to cancel out at short intervortex 
separations, whereas there will be interactions at large 
intervortex separations that are mediated by the gauge 
field. 

In the general model with the mixed gradient terms 
considered here (i.e. with pd ^ 0), there is in addition 
unscreened 1 /r interaction between vortices belonging to 
different condensates which is mediated by the direct 
Andreev-Bashkin drag. Thus, contrary to the Pd = 
case, there will be unscreened Coulumb interactions at 
all length scales. 



Observe that in the limit, ei = 62 = 0, Eq. 



(|18|) elimi- 
ting to 



nates Yukawa-type interaction potential and resu 
only only long-range interactions V{q^) = R/q^ like in a 
two-component superfluid with Andreev-Bashkin effect, 
see Ref. [13l Observe also that in contrast to the neu- 
tral model in Ref. [T31 in the above case when ei 2 7^ 
one always has a bound state of vortices which has finite 
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energy per unit length, as discussed in Sec. II B 

Thus, the vortex-vortex interaction matrix shows that 
adding the mixed gradient Andreev-Bashkin-type drag 
term to a superconductor, where components interact 
only via a fluctuating gauge field, might significantly al- 
ter the physics of fluctuations as a consequence of a sub- 



stantial change of the interactions between topological 
excitations. 



Finally, in the spirit of Section |II A[ we may rewrite 
the action in (17) in a form where the charged and the 



neutral modes are explicitly identified. 



S = 2Ptt' / d-^q 



P1P2 - Pd{pi + P2) 



[e2mi(q) - eim2(q)] ^ [e2mi(-q) - eim2(-q)] 



[ejRjimi{q) + eji?j2m2(q)] -1^ — 2 [efci?fcimi(-q) + efei?fe2m2(-q)] 



(20) 



Note that the vortex fields in the neutral sector inter- 
acts by an unscreened Coulomb interaction only, while 
the vortex fields in the charged sector interacts by a 
screened Coulomb (Yukawa) interaction. From this it 
follows that the corresponding propagators are given 
by ([e2mi(q) - eim2(q)] • [e2mi(-q) - eim2(-q)]) 

and ([eji?jimi(q) + eji?j2m2(q)] • [eji?jimi(-q) + 
eji?j2ni2(— q)]) ^ q^ -\- TOq. Here, toq is the effective dy- 
namically generated gauge mass that is nonzero in the 
low-temperature phase and vanishes at the charged criti- 
cal point. Moreover, there is also a neutral critical point 
associated with ordering the neutral sector of Eq. (201, 



with a corresponding non-analytic variation in the tem- 
perature dependence of the coefficient of the g^-term. 

Note that for any value of pd, the interactions of the 
vortex fields in the neutral sector are independent of any 
variation in the charges ei and 62 provided that the ratio 
62/61 is kept fixed, as readily seen by inspection of Eq. 
(20). On the other hand, the interactions in the charged 



sector depends on the value of the charge in the Yukawa 
factor l/(g^ + Wq)- 

Given the very different form of intervortex interac- 
tions produced by the gauge field coupling and by the 
Andreev-Bashkin drag, the interesting case when these 
interactions compete with each other cannot be mapped 
onto the previously studied regimes of systems interact- 
ing only by gauge field or only by intercomponent drag. 
Investigating the physics arising from this competition is 
the main objective in this paper. 



III. DETAILS OF THE MONTE CARLO 
SIMULATIONS 

Large-scale Monte Carlo simulations were performed in 
order to explore the phases and phase transitions of the 
model ([I]). We discretize space into a three-dimensional 
cubic lattice of size L x Lx L with lattice spacing a — 1. 
The phases are defined on the vertices of the lattice, 
0j{r) — 6'rj and the phase gradient is a finite dif- 



ference of the phase at two neighboring lattice points, 
d^9j{r) — > Af^Orj = Or+p,j — 6r,j- As for the finite dif- 
ference of the phases, the gauge field is also associated 
with the links between the lattice points, A^{r) — ?> A^^^. 
Moreover, the curl of the gauge field yields a plaque- 
tte sum (V X A(r))^ J2^n^t^'^v^'^^r,,T Here, e^^^ is 
the Levi-Civita symbol. The compact phases O^j have 
to be 27r-periodic. This is accommodated by the Vil- 
lain approximation of the effective Hamiltonian,!^ which 
also yields a faithful lattice representation of the direct- 
current current interaction (i.e. drag) term.^i^! Our effec- 
tive lattice model thus reads 



Z 



VA 



-/3ff[ei,e2,A;/3] 



where the local Villain action is 



^=2 



(22) 



Here, Ur,^,j ^fi^r.j ^^-^r,p 



- 2?™ 



r^fj,j IS a one compo- 



nent Villain argument. The sum over the integer valued 



fields, 



, is from ~oo to cxi ensures 2n periodicity 



of the Hamiltonian with respect to the gauge invariant 
phase difference. 

All Monte Carlo simulations start with an initializa- 
tion of the system, either disordered, when all phases 
and gauge fields are chosen at random, or ordered, when 
phases and gauge fields are chosen constant throughout 
the system. Subsequently, a sufficiently large number of 
sweeps is performed in order to thermalize the system. As 
a valuable check on the simulations, the calculated quan- 
tities should be invariant with respect to the initialization 
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procedure. A Monte Carlo sweep includes local updat- 
ing of all five fluctuating field variables (compact phases 
G [0, 2n) and the non-compact gauge field Aj.^^) at all 
lattice sites in the system, according to the Metropolis- 
Hastings algorithm. There is no gauge fixing involved, 
as summation over gauge equivalent configurations will 
cancel out when calculating thermal averages of gauge 
invariant quantities. Moreover, periodic boundary con- 
ditions are applied in all simulations. 

In most cases, we also apply the so-called parallel tem- 
pering algorithm,'221 allowing a global swap of configura- 
tions between neighboring couplings, after the local up- 
dating is finished. The explicit temperature dependence 
in the Hamiltonian of the Villain modeP^ must be con- 
sidered when calculating the probability of exchanging 
configurations between two coupling values /?, /3', which 



We expect no extra singular behavior due to the tem- 
perature dependence in the Villain Hamiltonian, so the 



WpT 



I, ifA<0, 
e-'^, ifA>0, 



(23) 



where A = p'{H[X-p'] - H[X'-p']) - p(H[X;li] - 
H[X'; /?]), and X, X' are the configurations at /3, /?' ini- 
tially. To increase the performance of the parallel tem- 
pering algorithm, the set of coupling values was selected 
according to the initialization procedure in Ref. |32, to 
yield approximately the same acceptance rate for the par- 
allel tempering move throughout the entire range of cou- 
pling values in the simulation. By introducing the par- 
allel tempering algorithm, the quality of the statistical 
output was substantially improved by reducing the auto- 
correlation time at critical points by 1-2 orders of mag- 
nitude compared with conventional Monte Carlo simu- 
lations with local updates only. Even in regions of the 
phase diagram where coupling intervals were too large 
for configurations to access all coupling values within a 
reasonable amount of MC sweeps, which is required to 
take full advantage of the parallel tempering method,'22l 
an improvement of the statistical output was achieved. 



A. Specific heat 

We measure the specific heat per site C^, by the energy 
fluctuations. 



(24) 



where the brackets denote thermal average with respect 
to the partition function in (211). In fact. 



this ex- 
pression is not quite right for the Villain model be- 
cause of the explicit temperature dependence in the 
Hamiltonian.'^ Generally, the specific heat is given by 
L^Cv = —0^dU/{dji), where the internal energy is given 
by [/ = ~d\x\ Z / {dj3).^^ Thus, the specific heat is written 

c.L^ _l (dm) _/dm)\y _&\m\ ,,,, 
op \ dp I) dp^ /• ^'^^ 



singular behavior in (251 should also be captured in the 
energy fluctuations of (24). Thus, we expect Eq. (24) 



to reproduce the correct critical behavior of the heat ca- 
pacity, as was the case in Ref. 34. In practice, both 
equations were used, and the results were identical with 
respect to critical behavior. In the analysis of the Monte 
Carlo simulations, the critical temperature of the phase 
transitions was determined by locating the anomaly of 
the heat capacity, and the same critical temperature was 
found with both equations. 



B. Helicity modulus 

The helicity modulus is a global measure of phase co- 
herence in a superfluid (i.e. decoupled from gauge field) 
order parameter. It measures the energy cost associated 
with an infinitesimal twist 8 in the phase of an order pa- 
rameter across the system. In order to obtain the correct 
energy cost with respect to composite phase combina- 
tions such as e.g. phase difference, we must perform a 
general twist in a linear combination of the order param- 
eter phases, 

Or. J d'r,j = Ov,j -ajS- r, (26) 

where aj now is a real number associated with the phase 
twist in component j. By selecting ai, 02, we may 
measure the phase coherence of any linear combination, 
fli^i +026*2, in order parameter space. That is, if we want 
to measure the helicity modulus of the neutral mode as- 
sociated with the phase difference we must impose a twist 
in the phase difference, i.e. Oi = 1, 02 = — 1. In general, 
the helicity modulus is given by the second derivative of 
the free energy with respect to the infinitesimal twist 



^ /^,(ai ,a2) 




(27) 



J (5=0 



where the notation 0' simply means that all phase vari- 
ables are replaced according to (26). In our case, with 
an isotropic system, we expect the helicity modulus to 
yield directionally independent results within statistical 
errors. For more details on the helicity modulus in the 
special case of the Villain model, we refer to Refs. Upland 

El 



C. Gauge mass 

To capture the properties of the gauge field A, we 
study the gauge field correlator (AqA_q), expHcitly 
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given for the lattice model 
2 



(28) 



where |QqP is the Fourier representation of the discrete 
Laplace operator, given by |QqP = sin(7r(7^/L)]^, 
and 



q,m/ 



(29) 



is the correlation function of the linear combination of 
vortex fields that corresponds to the charged sector of Eq. 



(20). Here mqj- is the lattice model vortex field of com- 



ponent j in Fourier space. The details of the derivation 
are given in Appendix [Bj In particular, we will use this 
quantity to extract the order parameter for the normal 
fiuid-superconductor phase transition, i.e. the dynami- 
cally generated gauge field mass, or Higgs-mass. The ef- 
fective gauge mass ttia is extracted from the gauge field 
correlator h^IMM 



q=0 /3(AqA_ 



(30) 



This quantity is employed as order parameter of the su- 
perconducting phase. Note that the dynamic creation 
of mass at Tc and the onset of the Meissner phase, the 
manifestation of the Higgs mechanism in London su- 
perconductors, is governed entirely by the long-distance 
behavior of the vortex correlator of the charged mode, 
cf. Eqs. (28) and (29). In the ordered phase, where 
vortex loops are confined, limq_j.o(mfc(q)m„(— q)) ~ 
g^, such that limq_j.o(A(q)A(— q)) ^ const, ren- 
dering the gauge-field massive. When vortex loops 
proliferate, limq_^o(mfe(q)mm(~q)) 7^ 0, such that 
limq_j.o(A(q)A(— q)) ^ l/g^, rendering the gauge-field 
massless. 

In the Monte Carlo simulations, the vortex fields of 
both species are extracted from the phase and gauge 
field distributions by considering the plaquette sum of 
the gauge invariant phase difference 



E 



(31) 



where the left hand side is the plaquette sum of the gauge 
invariant phase difference, A^^rj" — ej^r,/^ and Tir./ij is 
the real space vortex field. The gauge invariant phase dif- 
ference must be kept in the primary interval for each link 
in the plaquette sum in order to accommodate vortices 
in the lattice model. Now, by Fourier transformation of 
the vortex field, C/c,q is calculated, and to find the gauge 
mass, curvefitting of the quantity 2/(/3(AqA_q)) is per- 
formed for small q-values in order to extract the q 
limit. 



IV. MONTE CARLO RESULTS, pd = 0, ei = 62 = e 

Here we present the simulation results for the case dis- 
cussed in Sec. |IIB| In this section we consider in general 
unequal stiffnesses pi ^ p2 in the regime where pd = 0. 

Fig. [ijshows the simulation results varying the stiffness 
P2, when the other stiffness pi is set to unity. Results are 
obtained for 6 different values of the electric charge, and 
we focus on the regimes where there is a strong compe- 
tition between proliferating topological defects. We lo- 
cate the critical inverse temperature of the charged and 
the neutral critical point by locating the anomaly of the 
heat capacity associated with the phase transition. The 
charged critical point is associated with the point where 
the Meissner effect sets in, evident by onset of the effec- 
tive gauge mass toa, whereas the neutral critical point 
is associated with the onset of the order in the gauge 
invariant phase difference with a corresponding nonzero 
value of the associated helicity modulus (1 



A. Topological excitations 

Consider now the case when the neutral critical line is 
situated above the charged critical line, that is, when go- 
ing from phase I (t/(l) x [/(I) broken symmetry) to phase 
III (broken C/(l) charged symmetry) across the neutral 
phase transition line in Fig. [T] This phase transition is 
driven either by proliferation of (0, 1) or (1, 0) vortices. 
The composite (1, 1) vortices do not couple to the neu- 
tral sector of Eq. ([9| and can thus never be responsible 
for destroying the order in the neutral sector. The other 
composite topological excitation (1, —1) is, by inspection 
of (|9|, seen to have neither energetic nor entropic ad- 
vantage over individual vortices. Because the vortices 
(1, 0), (0, 1) cost the same amount of energy in the neu- 
tral sector, but the vortex with lowest stiffness pj costs 
less energy in the charged sector, the neutral critical line 
must be associated with proliferation of individual vor- 
tices of the component with the smallest value of the 
bare stiffness pj, when going from phase I to phase HI. 
This phase transition is, outside the region where there 
is a strong competition between different kinds of vortex 
excitations, found to be of second order in the 3Dxy uni- 
versality class.l^ When the individual vortices proliferate, 
the corresponding stiffness is renormalized to zero and 
the remaining condensate will be a charged condensate 
with order in the remaining component. Thus, the re- 
maining condensate will, at a higher temperature, have 
a phase transition similar to that of a one-component 
superconductor, 

Fni-.iv - |(Ve, - eAf + (32) 

where j now is the index of the component with largest 
stiffness pj. 

This is verified in Fig. [T] by observing that the 
charged critical line between HI and IV asymptotically 
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P2 P2 

Figure 1: (Color online) The phase diagram in the (p2, /3)-plane for the model in ([9| at 6 increasing values of the electric charge 
e when p\ = 1. Blue x markers (x) connected with dashed lines are charged critical points, and red squares (□) connected with 
dotted lines are the neutral critical points. When these critical points are merged, it is shown by filled squares (■) in violet 
connected with solid lines. Moreover, with lines in silver color, we present critical points of one-component superconductors 
with e as denoted in panel. The horizontal line is the critical line when p = p\ = 1, and the plus markers (+) are the critical 
points when p = p2. For these reference lines, the dashed and dotted line type correspond to charged and neutral critical 
points, as above. The inset in panel (b) is a magnification of the lines merger region. Phases are denoted by roman numbers. 
I. Ordered phase with spontaneously broken (7(1) x U{1) symmetry, mA 7^ 0,T^^(i__i) 7^ 0. II. Spontaneously broken global 
(7(1) symmetry, with restored (7(1) gauge symmetry, mA = 0,Tp,(i _i) 7^ 0. III. Spontaneously broken (7(1) gauge symmetry, 
with restored global (7(1) symmetry, mA 7^ 0, Tp.(i__i) = 0. IV. Normal phase with fully restored (7(1) x (7(1) symmetry, 
mA ~ 0, T^^(i,_i) — 0. The system size considered is 32"^. Except for inset, error bars are smaller than marker size and thus 
omitted from diagram. 
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approaches the one-component reference hnes away from 
the region of competition between different kinds of topo- 
logical excitations. This phase transition is second order 
of the inverted 3Dxy universality class. ^ 

When there is a transition from phase I {U(l) x U{1) 
broken symmetry) to phase II (broken U{1) neutral sym- 
metry) in panel (e) and (f ) of Fig. [ij the charged critical 
point is situated at a lower temperature than the neu- 
tral critical point. In this case the topological defects 
responsible for the phase transition are (1, 1) vortices, 
because the other possible vortices will destroy order in 
the neutral sector of Eq. and thus are not proliferat- 
ing at this transition line. As discussed in Sec. II B the 
(1,1) vortices proliferating from an ordered background 
may be mapped onto a single component superconductor 
with effective stiffness (if to neglect their internal struc- 
ture) p' = pi+ P2, 



7eff 
I^II 



Pi + P2 



{Ve - eA) 



(V X A) 



(33) 



In Fig. [2] we show results when bare component stiff- 
nesses pj are kept fixed, and electric charge e is varied. 
In panel (a), we also present a one-component reference 
line corresponding to the phase transition of the super- 



conductor in (33). Indeed, away from the splitting point. 



the transition from I to II approaches this reference line, 
and the actual second order phase transition is of the in- 
verted 3Dxy universality class. Note that the mapping 
in ( 33 ) yields a one-component superconductor with stiff- 
ness pi + p2 that always is stiffer than the two reference 
lines in Fig. [l] (which are one-component superconduc- 
tors with stiffness pi and ^2)- Thus, the charged tran- 
sition line between phase I and phase II is always lower 
than the reference lines in Fig. [ij Phase II in Fig. [ijand 
[2] is the metallic superfluid phase (i.e. exhibiting order 
only in the gauge invariant phase difference) discussed in 
Sec. |IIB and the effective free energy in the remaining 
superfluid condensate is given in Eq. (11). The cheap- 



est topological defects that proliferate at higher temper- 
atures and destroy the remaining composite order in this 
phase, are individual vortices. Hence, away from the re- 
gion of competing topological defects (i.e. away from the 
splitting point) , the transition line from phase II to phase 
IV should be similar to a one-component superfluid with 
effective stiffness p' = piP2/{pi + P2), 



F\ 



P1P2 



II^IV 



2{pi+P2) 



{vef 



(34) 



Note that in both panels of Fig. [2] the neutral tran- 
sition line between II and IV is found to be asymptoti- 
cally independent of e, thus approaching a constant value 
asymptotically far away from the region of competition 
with different vortices, as Eq. (34) suggests. Moreover, 
Eq. (34) predicts the value 15 = (ic{pi + P2)/{piP2) of 
the actual line, which corresponds well with the results 
in the figure. Here, /3c ~ 0.334 is the critical point of the 
one component superfiuid (e = 0) when p = 1. 



Note that vortices on the form (ni,n2) with Uj > 1, 
n-k^j > 1 can, by inspection of (|9]), be shown to always be 
energetically unfavourable compared with other topolog- 
ical excitations in this model. Such higher order vortices 
are thus not relevant when pd — and ei = 62 = e. 



B. Gauge field fluctuation driven merger of the 
phase transitions in case of unequal bare stiffnesses 



We next discuss the evolution of the phase diagrams 
m Fig. [l]and[2]whe n e is varied. When charge increases, 
the energy of the composite (1,1) vortices (which have 
no topological charge in the neutral sector), as well as 
the energy associated with charged currents of individual 
vortices decrease. This leads to a formation of a region in 
the phase diagram which is characterized by a merger of 
the two U{1) transitions in the case of unequal bare stiff- 
nesses of the two condensates. Thus, even in the case 
of unequal stiffnesses, when the coupling to a fluctuating 
noncompact gauge field is sufficiently strong, there ap- 
pears a phase transition directly from the ordered phase 
with spontaneously broken C/(l) x U{1) symmetry to the 
fully disordered normal phase. See also discussions of 
transition mergers caused by other kinds of couplings in 
Refs. 13 14 36, Panel (b) of Fig. [2] clearly illustrates 
this behavior. In this panel, the value of bare stiffness 
disparity is fixed when e increases. For low values of e 
there are two phase transitions: At lower temperature 
individual vortices with lower stiffness proliferate while 
at higher temperature a proliferation of individual vor- 
tices of stiffer condensate takes place. However, when e 
increases, the two lines approach each other and merge 
at e ss 1.3. 

The line merger is a consequence of the fact that at 
a substantially large electric charge, the bare energy of 
an individual vortex in a broken U{1) x U{1) phase is 
dominated by the neutral mode. Because a proliferation 
of less energetically expensive individual defects destroys 
the neutral mode, this eliminates the bare long-range log- 
arithmic interaction between vortices in the stiffer con- 
densate, leading to a dramatic decrease in their bare 
line tension and thus to their preemptive proliferation. 
On the other hand in a range of parameters a prolifera- 
tion of composite (1,1) vortices can trigger proliferation 
of individual vortices again leading to a "preemptive" 
restoration of the full U{1) x U{1) symmetry via a sin- 
gle phase transition. When electric charge is increased 
further, then eventually at a certain point in the inter- 
val e e (1.75. ..1.875) the (1,1) vortices become much 
less energetically expensive than other excitations and 
can proliferate at low temperatures without triggering 
a proliferation of individual vortices. Then the metallic 



superfluid phase (II) emerges as discussed in Sec. II B 
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e e 

Figure 2: (Color online) The phase diagram in the {(3, e)-plane for the two component 3D London model Q with pi — 1 
and for two different values of p2- In the left diagram is p2 = pi = p = 1, whereas p2 — 1.15 in the right diagram, i.e. there 
is a moderate disparity. Markers and line types are the same as in Fig. [l] i.e. blue x markers (x) connected with dashed 
lines are charged critical points, red squares (□) connected with dotted lines are neutral critical points, filled squares (■) in 
violet connected with solid lines are merged transitions. The silvered plus markers (+) in the left diagram is a one-component 
reference line of a superconductor with bare stiffness 2p and charge e. Roman numbers denote the different phases as given in 
the caption of Fig. [T] Note that these diagrams are 2D cross sections of a 3D phase diagram in p2, e)-space perpendicular 
to the cross sections in Fig. fl] The lattice size is 32^. Errors are smaller than marker size and thus omitted from diagram. 



Order of the phase transition associated with 
the merged lines 



Competing gauge field and drag interactions in 
the case pi = P2 = 1 



Let us now characterize the phase transition along the 
merged lines of Fig. [l] and [2j In Ref. [ini using the j- 
current model the transition line from C/(l) x J7(l) to fully 
symmetric state in the case of equal stiffnesses presented 
in panel (a) of Fig. |2] was found to be a first order 
transition. We obtain consistent results in our Villain- 
model based simulations. 

Furthermore in Fig. [3] we report the simulation re- 
sults associated with the merged line in a case when bare 
stiffnesses are not equal. We find a first order transition 
along the merged line in our case when there is a dis- 
parity of the bare phase stiffnesses. This shows that the 
first order phase transition ina[/(l)x[/(l) non-compact 
gauge theory is not related to the specific degeneracy of 
the model with equal stiffnesses pi = P2, but appears to 
be related to the case when there are several competing 
or composite topological defects. 



V. MONTE CARLO SIMULATION, GENERAL 
MODEL WITH BOTH GAUGE FIELD AND 
DISSIPATIONLESS DRAG INTERACTIONS 



Next, we present results from Monte Carlo simulations 
when both drag and gauge field mediated interactions are 
included. 



In Figure |4j we present results for the case when the 
bare component stiffnesses are equal pi = p2 = P = 1, 
and the gauge field couplings are equal, ei = 62 = e. We 
vary the inverse temperature l3 and the bare drag coeffi- 
cient pd and map out the phase diagram in the {P,l3pd)- 
plane for a number of different values of e. We consider 
positive Pd only. In this specific case, the charged and 
neutral modes in ([5]) are written, 



F = 



P - 'ipd 



(V6'i + V6I2 - 2eA)2 + (V X A)^ 



(35) 



Here, we have the interesting situation where drag- and 
gauge field mediated intercomponent long-range vortex 
interactions are found to be of opposite signs, see Eq. 
(18). Thus, the drag coupling pd 0, when significantly 
strong, favors formation of the (1, —1) composite vortices 
(via a mechanism similar to that in Ref. ITS). On the 
other hand, the gauge field coupling favors the formation 
of (1, 1) bound states of individual vortices when ei and 
62 are of the same sign. This competition is studied in 
Fig. |4] Its most striking consequence is that it leads to 
the existence of four phases: At strong drag there is a su- 
perconducting phase where a neutral mode is destroyed 
by the proliferated (1,-1) vortices (Phase V). At strong 
electric charge there is a superfluid phase with prolifer- 
ated (1,1) vortices (Phase II). 
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Figure 3: (Color online) (a): Histograms for the probability 
distribution of the internal energy per site U/L^ at /3 « 0.487 
when parameters are pi = 1, p2 = 1-05 and e = 1.5. This 
is the merged transition point found in Fig. [T] at (p2, /9) = 
(1.05, 0.487) in panel (c). A double-peak structure develops 
when L increases, (b): Upper panel shows that the finite size 
scaling of the latent heat per site A[/ /L^ approaches a finite 
value when L increases. This is the distance between the 
peaks in (a). The lower panel shows the finite size scaling of 
the difference in the free energy, A_F = (1//3) ln(Pniax/fmin) 
taken between the double peak value Pmax and the value of the 
minimum in between Pmin of the histograms in (a). For a first 
order phase transition, A_F ~ L''"^ .I^S Ferrenberg-Swendsen 
reweighting was used to obtain histograms with similar height 
peaksPl 



We next consider these phases more closely. The re- 
sults in Fig. |4] show that the phase V appears when 
(1,-1) vortices proliferate and thus there is no longer 
a broken symmetry in the neutral sector of Eq. ( 35 1 . 



Note that when we are well above the region of com- 
peting topological defects in Fig. |4] then, by neglecting 
the internal vortex structure, we may approximate the 
(1,-1) vortices to map onto vortices in a one component 
superfluid with stiffness p' = 2{p — 2pd), 



off 

I^V 



{p-2pa){vef 



(36) 
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Figure 4: (Color online) Phase diagram in the (/3, Ppd) plane 
with competing gauge field and drag interactions. Results are 
given for five different values of the charge e\ = e2 = e. The 
bare component stiffnesses are equal, pi = p2 = 1, and the 
system size considered is 32^. The gray-shaded area is the 
prohibited region by the stability condition, Eq. ([8|. Line 
type corresponds to character of phase transition as in Fig. 
[ij that is, charged lines are dashed, neutral lines are dotted 
and merged transition lines are solid. Except for this, lines 
are guide to the eye, only. The inset in the figure shows 
the schematic structure of the phases in the diagram for all 
cases with e > 0. Roman numbers denote phases. Phases 
I, II and IV are the same as given in the caption of Fig. [T] 
whereas V is, similar to phase III in Fig. [l] a phase with spon- 
taneously broken (7(1) gauge symmetry, and restored global 
17(1) symmetry, ttia / 0, T^_(i__i) — 0. However, in V the 
broken (7(1) gauge symmetry is associated with composite 
phase sum, whereas in phase HI of Fig. [l] it is associated 
with the phase of the single ordered component. For the given 
ranges of the phase diagram, II is only found for e = 2.25 
and e = 3 and V is not found for e = 3. When e = 0, all 
phase transitions are neutral and phase II and V are associ- 
ated with broken global ?7(1) symmetry in the phase differ- 
ence and phase sum, respectively.^ The results for e = are 
here simulated with a fluctuating gauge field, and coincide (as 
they should) with the equal stiffnesses results in Ref. 1 131 with 
no fluctuating gauge fleld. 



this physics manifests itself in the fact that in Fig. |4j the 
actual transition is seen to approach asymptotically the 
reference line Pp^ = {j3 — /3c/2)/2. 

In Sec. |IV[ the superconducting phase III, which sim- 
ilarly to phase V exhibits charged order and neutral dis- 
order, was created from the fully ordered phase by prolif- 
eration of individual vortices when we increased dispar- 
ity in the bare stiffness of the two components. Here, 
phase V is created by proliferation of composite vortices, 
and the coupling constant responsible for creating the 
phase is pd- Consequently, the remaining order is now 
in the gauge invariant phase difference of the charged 
mode, given by second and third terms in ( 35 ) . On the 



This effective limiting model is e-independent. Indeed 



other hand, phase III exhibits order in the gauge invari- 
ant phase difference of the component with largest bare 
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stiffness. Note also tfiat in tfie U{1) x (7(1) (Phase I) state 
with equal stiffnesses the (1,0), (0, 1) vortices carry half 
of the magnetic flux quanta-l^Sl It can be seen from ([35]) 
that in the phase V (1,0), (0, 1) vortices become equiva- 
lent and no longer have logarithmic divergence of internal 
energy per unit length due to absence of a neutral mode. 
I.e. they become similar to Abrikosov vortices, but carry 
only a half quantum of magnetic flux. This phenomenon 
is related to the fractionalization of superfluid velocity 
quantum in the metallic superfluid state.^^ From ( 35 ) it 



also follows that the individual vortices behave as vor- 
tices in a one component superconductor with effective 
stiffness p' — p/2 and double effective charge e' = 2e, 



v^iv 



P 



(ye - 2eAf 



(V X A)2 
2 ' 



(37) 



In Fig. [4] the transition from the phase V to the normal 
phase IV is indeed found to tend asymptotically to a 



phase transition one would predict from the model ( 37 1 
For this model, the transition line is found to be vertical 



in accordance with the drag independent stiffness in (37 1 . 
Note that when e increases, the critical temperature of 
the vortex loop proliferation is decreased and the vertical 
line moves to the right in Fig. |4j 

Next, the phase II may be investigated in a similar 
way as the phase V above. Phase II appears when (1, 1) 
vortices proliferate. As discussed in section |IV[ the re- 
maining order is in the neutral sector of (I5| and the 
transition to the normal state is governed by prolifer- 
ation of individual vortices that asymptotically behave 
as a one-component superfluid with effective stiffness 
p' = (p-2prf)/2, 



^II^IV 



P- 2pd 



(38) 



The critical phase transition of this condensate will follow 
the line jSp^ — {/3 — 2/3c)/2. Indeed, this is the case for 
e = 3 in Fig. [4] away from the region with competing 
topological defects. 

Similarly to Sec. IV we flnd evidence of a first or- 
der transition when lines are merged and e > 0, as seen 
in Fig. [5] When only drag or gauge field is included 
in a two-component system, first-order transitions may 
emerge.'^Sinilll Q^r results show that the first-order char- 
acter of this phase-transition line persists also in the case 
where both of the interactions are present and compet- 
ing. 



B. The regime where gauge field and drag 
interactions both favor formation of similar paired 
phase 

In Fig. [6] we present the phase diagram in the case 
when pi = /C2 = P = 1 and ei = —62 = e = 1. The 
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Figure 5: (Color online) (a): Histograms for the probability 
distribution of the internal energy per site U/ at the critical 
point when parameters are pi = p2 = p = 1, ei = 62 = e = 1 
and pd = 0.4. This is a merged transition point at (/3, Ppd) 
= (0.948, 0.379) along the critical line for e = 1 in Fig. [t] A 
pronounced double-peak structure is found to develop when L 
increases, (b): Upper panel shows the finite size scaling of the 
latent heat per site AU/L^. This is the distance between the 
peaks in (a). The lower panel shows the finite size scaling of 
the difference in the free energy, AF — (1//3) ln(Pniax/f'min) 
taken between the double peak value Pmax and the value of the 
minimum in between Pmin of the histograms in (a). For a first 
order phase transition, A_F ~ L'*"^.'^ Ferrenberg-Swendsen 
reweighting was used to obtain histograms with peaks of sim- 
ilar heightpl 



separation in neutral and charged modes is now, 
P - 2pd , 



F 



-{VOi - V6I2 - 2eA)^ + (V X A) 



(V01+V02)' 



(39) 



The motivation for investigating this particular case is 
found in the off-diagonal elements of the matrix in Eq. 
(18 1 where the interactions originating with the gauge 



field will act in unison with the bare drag interactions 
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upon switching the sign of the electric charge in one of 
the components (in contrast to the situation considered 
in the previous section). Indeed such a regime should 
lead to the growth of the paired (MSF) phase in the phase 
diagram. 

Consider the simulation results shown in Fig. [6j For 
comparison, we include the results when there is no gauge 
field coupling, ei = 62 = 0, and when gauge field coupling 
competes with the drag interaction, ei = 62 = 1. First 
notice that the paired phase which appears when charges 
are opposite, is the metallic superfluid phase (VI) which 
now is associated with spontaneously broken global U{1) 
symmetry in the phase sum (and not the phase difference 
as in Fig. [T] [2]and|4]). Positive drag will favor (1,-1) 
vortices as before. However, because of the change of 
sign of one of the charges, the (1,-1) vortices are now 
associated with the charged sector of Eq. ( 39 1 . The 



(1, 1) vortices are associated with the neutral mode, and 
thus the neutral critical point is determined by the onset 
of the associated helicity modulus (1^1). The gauge 
field renders the (1,-1) vortices the topological objects 
with lowest excitation energy. When they proliferate the 
superconducting sector is destroyed. Asymptotically, the 
associated phase-transition line is therefore expected to 
follow the behavior of a one-component superconductor 
with p' = 2{p — 2pd) and effective charge e, 



2pd){Ve-eAf 



(V X Af 



(40) 



The remaining condensate will have superfluidity de- 
stroyed via proliferation of individual vortices which 
asymptotically can be mapped onto a one-component su- 
perfluid with stiffness p' — p/2^ 



(41) 



his is the exact same behavior as expected when e = 0, 
which is also confirmed by simulations in Fig. |6] Note 
that there is neither pd nor e dependence of this line. 

Fig. [6] also shows that when gauge field and drag 
act in unison it amounts to a small increase of the re- 
gion of paired phase compared to the case when there is 
only drag interaction. However when interactions com- 
pete there is a stronger effect of the suppression of the 
corresponding paired phase. Also note that the cases 
ei — —62 = 1 and ei = 62 = 1 coincide when pd = 0. 
This is readily inferred from Eq. [T] since when pd — 0, 
one model can be mapped onto another by change in sign 
of charges accompanied by a sign-change of one of the 
phases V6j. 



VI. CONCLUSION 

We have considered a three-dimensional lattice super- 
conductor model in the London limit, with two individ- 
ually conserved condensates. These condensates interact 
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Figure 6: Phase diagram in the (/?, fipd) plane of the general 
model when p\ = p2 = ^ for the case of different charges 
ei = —62 = 1. These are the black x markers (x) and the 
line type denotes charged, neutral and merged critical lines by 
dashed, dotted and solid lines as in Fig. [T] [2] and |4] Roman 
numbers denote the phases of this particular case, ei = —62 = 
1. I. Ordered phase with spontaneously broken (7(1) x ?7(1) 
symmetry, mA 7^ 0, T^_(i,i) 7^ 0. IV. Normal phase with fully 
restored C/(l) x (7(1) symmetry, mA = 0,Tp_(i i) — 0. VI. 
Spontaneously broken global [/(I) symmetry, with restored 
?7(1) gauge symmetry, rriA = 0, T^.(i_i) 7^ 0. For comparison, 
the results of the two cases ei = 62 = 0, 1, from Fig. |4] 
are presented. The phases for these two cases follow from 
the inset and caption of Fig. |4] The lines are guide to the 
eye. The system size considered is 32^. The uncertainties in 
the position of the phase-transition lines are smaller than the 
marker size, and are omitted from the diagram. 



with each other by two mechanisms. The first is a dis- 
sipationless Andreev-Bashkin drag term representing a 
current-current interaction. The second is a fluctuating 
gauge field. Intercomponent Josephson coupling is ab- 
sent on symmetry grounds. Such models are relevant 
in a number of physical circumstances ranging from the 
theories of the quantum ordered states of metallic hydro- 
gen, models of neutron stars, and were earlier suggested 
as effective models describing valence bond solid to Neel 
quantum phase transition in the proposed theories of de- 
confined quantum criticality. 

In the case when there is no intercomponent drag, pd — 
0, and component charges are equal, ei = 62 = e and 
there is a disparity of the bare component stiffnesses, we 
find that a sufficiently strong coupling to a non-compact 
gauge field causes a merger of phase-transition lines. This 
yields a direct transition from broken [/(I) x [/(I) to nor- 
mal state even when the bare component stiffnesses are 
unequal. When the charge e is increased, the merger 
occurs for a higher disparity of stiffnesses. However, a 
further increase of the coupling beyond a certain critical 
strength results in a new splitting of the transition line. 
This yields a metallic superfluid phase. The merger of 
the [/(I) transition lines is associated with a competition 
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between different kinds of topological defects where pro- 
liferation of one type of vortices triggers a "preemptive" 
proliferation of another. The result is a much more com- 
plex picture of the behavior of topological defects in the 
phase transition than in singel-componcnt U{1) models. 
The second splitting is due to the fact that increased cou- 
pling to the non-compact gauge field decreases the free 
energy per unit length of a bound state of topological de- 
fects. The bound state in question (a composite vortex) 
has a topological charge only in the charged sector of the 
model. This in turn results in increased suppression of 
the critical stiffness associated with the charged sector 
of the theory, which eventually undergoes a symmetry- 
restoring phase transition before the neutral sector. 

We find that also when the bare stiffnesses are un- 
equal, the merged phase transition is first-order in char- 
acter. Note that previously first order transitions were 
reported in t he U{1 ) x U{1) gauge theory with degener- 
ate stiffnesses j^" ! ^"'^ ! U{1) x U{1) models with a compact 
gauge-field, as well as to phase transitions in the SU (2) 
model with noncompact Abelian gauge field.'i^ 

At stronger coupling, far away from the region of 
merged lines phase diagram exhibits 3Dxy (neutral) or 
inverted 3Dxy (charged) critical points. 

In the main part of the paper, we have performed a 
study of the phase diagram of the generic C/(l) x [/(I) 
lattice London gauge model featuring both gauge field 
and direct non-dissipative drag interactions. We have 
obtained, through large-scale Monte Carlo simulations, 
its phase diagram as a function of these two generic cou- 
pling constants. 

For the case where the bare component stiffnesses and 
charges are equal, pi = P2 — ^ and ei — 62 = e, we 
find the formation of two different paired phases as a re- 
sult of a competition between gauge-field and intercom- 
ponent drag couplings. High values of drag produce a 
composite superconducting phase associated with a bro- 
ken local U{1) gauge symmetry in the phase sum. There, 
the theory effectively features a doubled electric charge 
compared with U{1) x U{1) phase, cf. Eq. (37). At high 



values of e, the gauge field coupling wins over the drag 
coupling, yielding a paired superfluid phase (the metallic 
superfluid) associated with the order in the gauge in- 
variant phase difference. In between these two different 
phases, there is a region with a direct transition from bro- 
ken U{1) X U{1) to normal state which exhibits clear-cut 
signatures of a first-order transition, cf. the transition 
line connecting regions II and V in Fig. [4j 

For comparison, we also reported a quantitative study 
of the situation where gauge-field mediated intercompo- 
nent interactions and intercomponent drag both favor 
metallic superfluid phase. 
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Appendix A: The A'^-component model 

We may easily generalize to the case of arbitrary num- 
ber of components N. Again, the action will be given on 
the form 



J VOi... J V6n J VA e"", 
S^^J rfV{[V X A(r)]2 

+ [V0,(r) - ejA(r)]i?,fc[V0fc(r) - efeA(r)]}. 
The matrix Rjk is in general given by 



(Al) 



Rjk - 



(^Pj - X! P'^->'^ + Pd,3k, (A2) 



where Pd.jk is the drag-coefficient between components 
j and k, obviously, = Pd,k] and pd.jk = when 

j = k. Following exactly the same procedure as in the 
case N = 2, we arrive at the A^-component action 



j ( - ^2 ^■'^2 ) Ufe(-q) 



■A(q)(g2+rn2)A(-q) 



(A3) 



where the Fourier transform of V0j(r), is denoted by 
Uj(q), and 



ejRjkek- 



(A4) 



This expression is seen to reproduce the case N — 2 
given in Eq. ([t]). The gauge field is integrated out and 
the dualization now follows the same path as previously, 
yielding 



Z 



mi niN 

S = 2^7r^ J Am, (q)F,fc(q2)mfe(-q) 
where the vortex interactions are given by 



(A5) 



1 



Rjk 



ijk 



Ij^tnk 



eiEmRljlimk 



(A6) 
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This is seen to be on precisely the same form as Eq. ( 18 ) 
for the case N = 2. Following Appendix B in Ref. |9J a 
dualization of the corresponding two-component lattice 
model in Eq. (21) may be performed to yield the exact 



same result as in ( A5 1 and ( A6 ) where the vortex fields 



now are defined on the vertices of the Fourier space dual 



lattice and q 



IQc 



Appendix B: Gauge field correlator 

By adding source term and Fourier transformation of 
the model in Eq. (Al), the generating functional for 



deriving the gauge field correlator reads 



Zj = J V6i... j V0N J VA e-", 

^2 J ^'H^ [J(q)A(-q) + A(q)J(-q)] 

+ [U,(q) - e,A(q)]i?jfc[Ufc(-q) - efcA(-q)] 

(Bl) 

where J(q) are the electric currents that couples linearly 
to the gauge field in the source terms. Sum over repeated 
indices is assumed. We now proceed similar to Sec. |II C| 
by completing the squares of the gauge field and integrate 
out the shifted gauge field A(q) = A(q) + (J(q)//3 — 
eji?jfeUfe(q))/(g^ +ml) which yields 

Zj = J VOi... j V0N e-^, 
13 f i J(q)J(-q) 



(B2) 



J(q)eji?jfcUfc(-q) + eji?jfeUfe(q)J(-q) 



U,(q) i?, 



+ ml) 
q^ + m\ 



Ufe(-q) 



We now employ the constraint V • J(r) — 0, i.e. the 
electrical currents are divergence-free, such that compo- 
nents parallel to q are unphysical. Thus, the physical 



components of J(q) in the first term of (B2) are pro- 



jected out with the transverse projection operator 



(B3) 



As discussed in Sec. II C we disregard the longitudi- 



nal part of Uj{q) and introduce the Fourier transformed 
vortex fields by (16 1. Thus, the generating functional is 
written 



mi m TV 

5*0 = 2Pn^ j d\ mj(q)-l/jfe(g2)mfc(-q), 
Si^ I £q< 



q'^iq^+ml) 
X [m,,j(q)J^(-q) - J^(q)m,,j(-q)] (B4) 

2P{q^+ml) [' 



where Vjk{q^) is given by ( |A6[ ) and e^nri is the Levi-Civita 
symbol. 

The gauge field correlators are derived the standard 
way by functional derivation of the currents 



(A^(q)A,(-q)) = 



1 



Zo5J^{-ci)5J,{ci) 



J=0 



5J^(-q)5Jy(q) 



(B5) 



j=o/ 



where Z^ = Zj\j^o = Y^mi ■■■ Emjv ^° ^^^^ brack- 
ets denote thermal average with respect to Zq. The func- 
tional derivation is easily performed by expanding the 
exponential in series and keep terms of 0{J^), the only 
terms that survives both derivation and J = 0, to yield 



(5J^(-q)(5J^(q) 



j=o 



5j €[ J^j k-^hn^ fiaP^i^jK 

X gQg^TO^,fc(q)m„,„(-q) (^6) 
I3{q^+ml)' 

The product s^^ap^ u-yK is evaluated by the determinant 



^av ^aj ^OLK 
S/Su Sp^ 5(3^ 



(B7) 



to yield 

(^^(q)A.(-q)) = 



Air'^ejeiRjkRim 



2\2 



/3(g2 + m^) q'^{q'^+ml) 
X (PT,M^mfe(q)m„(-q) 
- »77,i.,fc(q)m^^„j(-q)), 



(B8) 



when (B6) is inserted in (B5). We now find the gauge 



field propagator by letting — > /i in ( B8 ) and summing 
over repeated indices, thus 



(A(q)A(-q)) 



47r2eje;i?jfei?/„(mfe(q)m™(-q)) 



(B9) 



+ 



Piq^+ml) 



2\ ■ 
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The gauge field correlator of the two-component discrete 
model in Eq. (21) is found similarly to Appendix C in 



Ref. ini and the result is as given in ( B9 ) with — > | Q 



and vortex fields defined on the vertices of the Fourier 
space dual lattice. 
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